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Abstract 

GRACE/SUSY is a program package for generating the tree-level amplitude and eval- 
uating the corresponding cross section of processes of the minimal supersymmetric 
extension of the standard model (MSSM). The Higgs potential adopted in the sys- 
tem, however, is assumed to have a more general form indicated by the two-Higgs- 
doublet model. This system is an extension of GRACE for the standard model(SM) 
of the electroweak and strong interactions. For a given MSSM process the Feynman 
graphs and amplitudes at tree-level are automatically created. The Monte-Carlo 
phase space integration by means of BASES gives the total and differential cross 
sections. When combined with SPRING, an event generator, the program package 
provides us with the simulation of the SUSY particle productions. 



1 Introduction 



It is widely expected that the supersymmetric (SUSY) extension of the stan- 
dard model(SM) of particle physics [1] is one of the most promising approaches 
that could resolve the problems remaining in the standard model. In partic- 
ular, the minimal supersymmetric extension of the standard model (MSSM) 
has been extensively studied in the last decade, because it has the simplest 
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structure containing the least number of particles, yet it is complex enough to 
allow the most essential features characteristic in any theory of SUSY. 

Over these ten years we have been working with the automatic computation 
of high energy physics processes. The system for the computation of the SM, 
GRACE, has been published in [2] . In order to be able to handle SUSY processes, 
we have made several necessary modifications and extension of the system 
[3,4]. First we have installed the new interaction vertices according to the 
MSSM lagrangian allowing an extension to a more general effective Higgs 
potential (section 2.4). Second we have modified the program so that we can 
deal with the fermion number non-conserving vertices and Majorana particles, 
both of which are inherent to any SUSY theory. 

Using this extended version of GRACE/SUSY, we have published a package 
of event-generator, SUSY23, which contains 23 specific SUSY processes for 
e + e~ — > 2-body and 3-body [5]. 

Since then we have made extensive checks of the system. In this paper we 
present the detailed description of the program system which enables us to 
compute automatically the cross section of any MSSM tree process, and ex- 
plain the structure of GRACE/SUSY together with instructions for how to use 
it. 

Program pacakge GRACE/SUSY (GRACE v2 . 2 . 0) is available form 
http : / /minami-home . kek . jp/. 



2 Theoretical framework of the MSSM 

2.1 Basic lagrangian 

The basic lagrangian of the MSSM in its full form has been given in sev- 
eral papers[6-8]. In this work we use the convention adopted in ref. [8], where 
the positively (negatively) charged charginos are considered as fermion (anti- 
fermion). The full form of the MSSM lagrangian is too lengthy to be repro- 
duced here. Readers should refer to the original work [8]. We show only the 
soft SUSY-breaking terms in order to define the SUSY breaking parameters. 
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-m^Hx - rh 2 2 W 2 U 2 - (m? 2 Hi ( ^ J ] H 2 + /i.e.) 
V 2m ^A tt H 2 A(g L )A( Mi? ) + ^MA d H 1 A(g L )A(rf R ) + h.c. 



V2 Vi 

+ y^h Ae H 1 A(£ L )A(e R ) + h.c 

-Y.[rfjlh + rn)Jy R ]. (2.1) 
/ 

In the first line $Hi and $h 2 stand for the iso-doublet Higgs superfields with 
hypercharge Y = — 1 and Y — +1 respectively, and in the second line A, A" 
and g a are the gaugino fields for U(l), SU(2) and SU(3) gauge symmetries, 
respectively. Hj appearing in the third to fifth line are the Higgs doublets 
which are decomposed as 



H2 " [%) = [{v 2 + ( t ) lUxl)/V2)- (2 ' 2) 

The scalar fields A(g^), A(u R ), etc. represent the sfermion fields, the super- 
partners of the corresponding fermion fields | ^ L j and u R , respectively, 

A(q L ) =[f L ), A(u R ) = u* R , A(d R ) = d R . (2.3) 



Although the so-called i?-parity violating interactions [9] are not included in 
this lagrangian and hence in GRACE/SUSY, users can implement them in the 
system by coding the necessary vertex functions etc. in appropriate subpro- 
grams. 



2.2 Particle content 



The MSSM consists of superfields containing the known fermions or bosons as 
their component fields, and two Higgs superfields. Note that in contrast to the 
SM, the MSSM includes three neutral Higgs particles and a pair of charged 
Higgs particles. The mass eigenstates of the CP-even neutral Higgs states are 
defined by 

0? 



[ H I _ f cosa sinaW^n , . 

U° J ~ I -sine* cosa J U° J ' 
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where the rotation angle a is determined in such a way that the Higgs mass 
matrix becomes diagonal. Since any of the known particles has its super- 
partner, the number of the elementary particles in the MSSM is more than 
double the number of particles appearing in the SM. In fact, for each known 
fermion / (except for neutrinos), there are two kinds of sfermions, the left- 
handed sfermion Jl and the right-handed sfermion The mass eigenstates 
of sfermion f\ with mass m; ( < m^) is defined by the orthogonal matrix 



f 2 ) [-sm6 f cos6>/ J [f R 



f = e,u,d,--- (2.5) 



where the rotation angle is determined by the mass diagonalization condition. 

For each gauge boson there is a gaugino, and to each Higgs particle corresponds 
a higgsino. Two neutral gauginos, A and A , and two neutral higgsinos, H® 
and L are combined to form four Majorana fermions called neutralinos, 

X® = | ^ j • The two-component spinors = 1 ~ 4) are linear combinations 
of A, A , Hi and H%, 

~\ = Vi(0 N ) ijXj , 2 = 1-4 (2.6) 

where Xj = H®, H®), and masses are ordered as m^o < m^o < m^o < 

m^o. The sign factor rji takes the value of +1 or +i, and is the orthogonal 
matrix which diagonalizes the mass matrix of the (A, A , H®, H®) system. 

The other neutral gauginos of the strong interaction, gluinos, remain unmixed 
with other particles. 

Corresponding to the charged Higgs particles if ± and charged gauge bosons 
W ± , there are two charged higgsinos, iff and ii^ 1- , and a pair of c harg ed 

r \ 

gauginos (wino), A 1 * 1 . Two charged fermions called charginos, xt = , 4 + 

I \l 

are linear combination of these particles, 



J [ - sin <f) R cos (j) R J [ Hi J ' 

A M = f 1 1 f cos< ^ sin ^l f A+ 1 (9 7) 

A^J l0 ej l-sin0 L cos <fi L J I Jf+ J ' 

The mixing angle <f)R, 4>l arid the sign factor are determined so that their 
masses satisfy the ordering < m-± < m-± Q. 

1 Note that the value of the mixing angles 4>L an d <pR become unstable when tan j3 = 
—{1/M2 or cot (3 = —fi/M2- If this happens, tan^x or tan^ become ±00. 
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Together with the particles appearing in the SM, these SUSY-particles con- 
stitute the ingredient of the MSSM. The list of the particles and their coding 
names is given in Appendix A. 

2. 3 Basic parameters of the MSSM 

The basic parameters of the MSSM come in when the lagrangian is expressed 
in terms of superfields [8] . They are 



coupling constants 


g,g',g s ,v 


gaugino mass parameters 


M 1 ,M 2 ,M 3 


fermion masses (for each generation) 


m e ,m u ,m d 


Higgs mass parameters 


m 1 ,rh 2 ,m 1 2 


trilinear couplings(for each generation) 




sfermion mass terms(for each generation) 


m e.Li m eRi m ULi m URi m d R 



The vacuum expectation values v i and v 2 are determined at the minimum of 
the Higgs potential and they are expressed in terms of these basic parame- 
ters. Assuming massless neutrinos, there are in all 4 coupling constants, 3Nq 
fermion masses, 6 + 8Nq soft supersymmetry breaking parameters, where Nq 
stands for the number of fermion generations. 



2.4 The MSSM parameters in GRACE/SUSY 



As alternatives to ten MSSM parameters for the coupling constants and gaug- 
ino and Higgs masses shown above, we can use in GRACE/SUSY the following 
set of ten parameters as inputs, 

a e , a s , M w , M z , tan/5, M A o,fi, M u M 2 , M 3 . (2.8) 



Unless specified explicitly by the user (see Appendix B), M\ and M 3 are 
computed from M 2 by the following unification conditions 



Mi = -tan 2 9 W M 2 , 
a 2 

M 3 = %sm 2 9 W M 2 . 



(2.9) 
(2.10) 
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Concerning the eight sfermion mass parameters in each generation, we use 
three mixing angles, and two up-type squark masses, the lighter down-type 
squark mass, and two down-type slepton masses as independent input param- 
eters. Due to the SU(2) L invariance of the £ so ft, the heavier down-type squark 
mass and the sneutrino mass are given as the solution of the following two 
equations (shown here only for the third generation sfermions), 



sin 2 9 b m~ = cos 2 9 t m\ + sin 2 9 t m~ 2 - cos 2 9 h m\ + m\ - m\ - cos 2(3, 

(2.H) 

ml T = cos 2 9 T m\ + sin 2 9 T m\ - m 2 + cos 2(3. (2. 12) 

Since we restrict our consideration only to the CP conserving processes, all 
the parameters are real. In particular, all the physical masses used as the 
input are real and positive. The allowed range of the SUSY parameters used 
in GRACE/SUSY whose values are not yet fixed by experiments are as follows 2 



-oo < n < +oo, 0<tan/3, < Mi, M 2 , M 3 , M A o, 

0<m fi <m h , -|<^<+|. (2.13) 

The CP-even Higgs, charginos and neutralinos are numbered as 

M h o < M H o, m-± < m^±, m^o < m^o < m^o < m^o. (2.14) 



The basic parameters of the MSSM are expressed in terms of these GRACE/SUSY 
variables given in (2.8) as follows, 



Ml 



M\ 



Aj = sin Of cos9f 



m~ — 171; 

h h 



fi cot (3, f 
\x tan (3, f 

2 T\/r2 



U, C, t 

d, S, b, 6, fi, T ' 



m f 

fn 2 j L = cos 2 9 fin 2 ^ + sin 2 9 fin 2 ^ — mj — M'^cos2f3(Ty — Qfs'fy), 

f = v e , z/ M , v T , e, ii, t, u, d, c, s, t, b, 



2 In ref.[7] the value of /i is restricted in the range < fi, while — oo < tan (3 < +oo 
is allowed. Due to the invariance of the lagrangian with respect to an interchange 
(n,tanP) <-> (—//, — tan (3), the physics with [i = /j,q > and tan/3 = — tan/3o < 
in ref.[7] corresponds to the one with fi = —fi < and tan/3 = tan/3 > in 
GRACE/SUSY. 
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m 2 j R = sin 2 Ofm 2 ^ + cos 2 dfin 2 ^ — in 2 — M§ cos 2f3Qf.s 2 v , 

f — e, /i,T,u,d,c,s,t,b. (2-15) 

In the last two equations of (2.15) the identities rh 2 L = m 2 e , m 2 L = fh\ , etc. 
are guaranteed by the relations (2.11) and (2.12). The vacuum expectation 
values V\ and v 2 are given by 



- v z w i — ■> 

z eyl+tan 2 /3 

> M w rvr 2 — tan ^ 



v 2 = 2-^ M 2 - M 2 '== . (2.16) 

M z eyl+tan 2 /? 

Note that the masses of Higgs particles, neutralinos and charginos are de- 
termined by the mass diagonalization matrices introduced in (2.4), (2.6) and 
(2.7), respectively. In particular, the Higgs masses are given at tree level by 



^HO,h° — 2 L 
M 2 



r 2 ±=M 2 + M, : 



'M% +M 2 z ± ^(M\ + M|) 2 - 4M 2 M 2 cos 2 2(3 
% + M 2 W . 



(2.17) 



According to (2.17), the MSSM at Born order predicts the existence of a CP- 
even Higgs particle which is lighter than Z. This possibility having already 
been experimentally excluded, we have to make a minimal modification of the 
Higgs potential so that M#o, M h o etc. can be treated as free parameters. For 
this purpose, inspired by the two-Higgs-doublet model [10,11], we introduce 
the effective Higgs potential (see also [12,13]) in the following way: 



V cS = CiH^Hi + C 2 U* 2 U 2 + 2C 3 Re(H 1 H 2 ) + 2C 4 [Im(H 1 H 2 )] 2 

+C 5 (H*H! - H;H 2 ) 2 + C 6 |HJH 2 | 2 + CMHtHx) 2 + C 8 (H;H 2 ) 2 

+4C 9 (HJH 1 )Re(H 1 H 2 ) + 4C 10 (H;H 2 )Re(H 1 H 2 ) 

+C n - (2.18) 

When the coefficients Cj take special values 



M 2 

C l =m 2 1 + /j 2 = M% sin 2 (3 cos 2(3, 



(2.19) 



M 2 



C 2 = m\ + /T = Mj cos^ (3 H cos 2/3, 

C 3 = m 2 2 = — sin f3 cos (3M\ , 
C,= l -{g 2 + g>\ 



(2.20) 
(2.21) 
(2.22) 
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C 6 = y, (2.23) 
C 4 = C 7 = C S = C 9 = C 10 = C u = 0, (2.24) 

the above effective Higgs potential reduces to the standard one of the MSSM. 
Since G\\ does not contribute to the Higgs masses nor to the interactions, the 
effective Higgs potential depends on ten parameters Ci, (i = 1 ~ 10), allowing 
more freedom than in the case of the original MSSM Higgs potential which 
contains three parameters, fh 2 , and rh\ 2 - The conditions at the potential 
minimum relate Ci to Vj, (j = 1,2). Therefore, the effective Higgs potential 
contains eight free input parameters 0, for which we take 

M A o, M H o, M h o, M H ±, tana, C 4 , C 9 , C 10 . (2.25) 



with the convention adopted in GRACE/SUSY 



~ < a < 0. (2.26) 



In terms of our input parameters, the coefficients Ci can be expressed as 



C 2 
C 5 

C 7 



cos 2a 



V2 



D sin 2a 



V2 
Vl 



C 3 , 



— [Mlo + ML - D cos 2a] — —D sin 2a - — C 3 , 
4 Av 2 v 2 

— sin (3 cos f3M\ + V\V%C± — vfCg — v lC w , 



D sin 2a 



3 M 2 A0 



4t> iv 2 2vf + v 2 



Cs 3 

r o°4, 



M'jr - Ml 



A 



vl + vl 



+ a 



vl 



vl 



V2 



Ci + -C 3 + (vf - v 2 2 )C 5 + 3v lV2 C 9 + 



Vl 

C 2 + ^C 3 

w 2 



''1 
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(vl 



vl)C 5 + -^C 9 + 3v 1 v 2 C 10 

V2 



where the vacuum expectation values Vi are given by (2.16) [f] and 



D = M% - M? 



h°- 



(2.27) 

(2.28) 
(2.29) 
(2.30) 

(2.31) 

(2.32) 

(2.33) 

(2.34) 



3 Compared with the MSSM, the effective Higgs potential (2.18) contains seven 
more parameters and consequently we need seven more new inputs. 

4 (2.32) and (2.33) are the conditions for the vanishing of the linear terms in the 
effective Higgs potential, thus assuring that v\ and V2 correspond to the minimum 
of the potential (2.18). 



8 



Table 1 

Input parameters used in GRACE/SUSY. 



coupling constants 




Gauge boson masses 


M z , M w 


Higgs sector 


M H o , M h0 , M A o,M H ± 
tan a 

C4, Cg, C10 


SUSY parameters 


H, tan/3, Mi, M 2 , M 3 


fermion masses 


m u , m c , m t 
m d , m s , m b 


sfermion masses 


ro di' m 5i' m W 


sfermion mixing angles 


#u> ^C) ^t) 

6<l, 6 S , Of,, 

6 e , 6^, 6 T 


particle widths 


Tw, Fz, ■■■ 



Combining (2.8) and (2.25), we use the following quantities as the inputs of 
GRACE/SUSY in the gauge and Higgs sectors, 

a e , a s , M w , M z , tana, tan (3, M A o, M H o, M h o, M H ±,fi, M 1 ,M 2 ,M 3 , 
C 4 ,C 9 ,C 10 . (2.35) 

In addition to these quantities, all the widths of the unstable particles must 
also be specified as input parameters. The set of the basic input parameters 
used in the program system are summarized in Table 1. The magnitudes of 
these inputs can be changed freely by the user. A list of the names of the 
variables created by the system in FORTRAN codes are given in Appendix B. 

A possible minimal extension of the MSSM is obtained by setting 

C 4 = C 9 = C10 = 0, (2.36) 

and using for tana and Mjj± the same form as given by the tree MSSM, 
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tan 2a = tan 2(3 



M 2 A0 + Mf 
M%>-M%' 



(2.37) 



M 2 H± = M 2 V + M 2 A 



(2.38) 



respectively. Concerning M H o and M h o, one can use any desired value, such 
as those evaluated by the one-loop MSSM. 

In GRACE/ SUSY, two modes are prepared for the choice of Higgs parameters. In 
mode 1 the strict MSSM Higgs potential is used, while in mode 2 the extended 
Higgs potential, (2.18), is used with the input parameters given in (2.25). The 
mode and the input parameters of the Higgs potential must be specified in 
subroutine sethptal. Note that in editing the subroutine, one can also use a 
different set of input parameters for the Higgs potential. See also section 4.4.2. 



2. 5 Treatment of the fermion number clashing vertices 

In models of SUSY the fermion number may not be conserved at the interac- 
tion vertices when they couple to neutralinos and charginos. Note that since 
the fermion number of chargino is not determined by their interactions, we 
define the positively charged charginos as Dirac fermions. 

In GRACE/SUSY we manipulate the fermion number clashing interactions by 
introducing the concept of G-line (GRACE-line), F-direction (fermion-direction) 
and M-direction (momentum-direction). The detail of the treatment for the 
tree case being given in ref. [14], we present here only a brief sketch of the 
procedure. 

Since the flow of a spin one-half particle is not disconnected nor has any 
branch, all the external spinors are connected pairwise. When such a spinor 
line connecting two external spinors is stretched in a straight line and displayed 
horizontally, we call it a G-line. According to the G-line rule to be discussed 
below, we compute the G-line amplitude corresponding to each G-line. The 
Feynman amplitude for each graph is then obtained by multiplying all the 
so-constructed G-line amplitudes, the bosonic parts which connect vertices in 
different G-lines and the relative phase of the graph at issue. 

G-line rule 

Let us assume that the process contains rij incoming spinors and rif outgoing 
spinors. We assign the particle number 1 to n« to the ni incoming spinor 
particles, and + 1 to Hi + nf to the rif outgoing spinor particles. For a given 
fermion line connecting two external spinor particles with particle number m 
and n ( m > n), we define a G-line in such a way that the particle m comes 
from the left and the particle n comes from the right endpoint. This defines 



10 



the G-line uniquely for each fermion line belonging to a Feynman graph and 
there are (rij + nf)j1 such G-lines at tree level. 



On each segment of the G-line separated by the interaction vertices, we assign 
two directions, an M-direction and an F-direction [15]. For the notational con- 
venience, we also introduce the notion of an A-direction, which is the direction 
in which spinors are ordered in the Feynman amplitude when a fermion comes 
in and goes out after some interactions, namely, the direction indicated by the 
left-oriented arrow. 

M-direction 

An M-direction is assigned in a natural way. For external incoming particles, 
their M-direction is the direction towards the vertex, while for the external 
outgoing particles it is the direction from the vertex. For the internal fermions, 
the M-direction is always taken in the A-direction (irrespective of whether the 
internal propagator is that of fermion or anti-fermion when seen in the A- 
direction). 

F-direction 

For Dirac particles we assign the F-direction in a natural way; for an external 
fermion the F-direction is the same as its M-direction, while for an external 
anti-fermion it is opposite to its M-direction. For the internal Dirac parti- 
cle it is taken to be equal to the A-direction for fermions (when seen in the 
A-direction) and to the opposite direction for anti- fermions. Concerning Majo- 
rana particles which are self-charge-conjugate, how to assign their F-direction 
is a matter of convention. In our convention, we define the F-direction of Ma- 
jorana particles by their M-direction. This means that the F-direction of the 
internal Majorana particle coincides with its A-direction and only the follow- 
ing type of the Majorana propagator appears in the G-line amplitude; 



where the momentum p runs from right to left. Clashing of the F-direction 
occurs only at vertices, and the charge conjugation matrix appears only at the 
fermion number clashing vertices. 

The above assignments of the F-direction and the M-direction uniquely iden- 
tify the external spinors as shown in Table 2. 

The general rule of flipping one of the F-direction at the vertex is summarized 
as follows 



1 




p 1 — to' 



Table 2 

Spinor assignment at the left and right endpoint of the G-lines. 



F-direction 


i — • 


— > • 


• i 


• > 


M-direction 












Hp) 


v{p) T 


u{p) 


v{p) T 




v(p) 


u{p) T 


v{p) 


u{p) T 



( • — )c~\ 

• — = ( • <— )C, 



where is a fermion line whose F-direction can be arbitrary. The charge 

conjugation matrix C satisfies the following properties 

& = C~\ C = -C T , C[1,75,7m,7^75] T C- 1 = [1,75,-7^7m7 5 ].(2.40) 

It is also important to note that we can choose the relative phase of the spinor 
wave functions such that 

Cv{k,X) T = u{k,X), Cu{k,Xf = v(k,X). (2.41) 

The rule of reversing the particle order at a vertex is given as follows; 

B — >• — ► A = -(A< — •< — B) T , 
B — >•< — A — —(A — >•< — B) T , 
B < — • — > A = -(A< — • — > B) T . 

Applying the above rules to the lagrangian of fermion number conserving Dirac 
fermions interactions (including chargino), 

C DD = ^(D a )T Da>D ^(Dp) + h.c, (2.42) 
one obtains the G-line rule for the Dirac particles (neglecting the overall factor 



D a < • < Dp r DaiD/3 , (2.43) 

Dp > • >■ D a -r T Da ,D - (2-44) 

Concerning the vertices including Majorana particles, corresponding to the 
lagrangian 
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C M d — * '(M a )T M a ,Dp^ '(D p) + ^(Dp)r DptMa 

*-MM = 

Q>/3 



(2.45) 
(2.46) 



where T^m^ = 7orj W/3i £ )ct 7o and Z? stands for a Dirac particle which can be a 



chargino, we obtain the following G-line rule, 



D a 

Mr 



D 



D 



Mp 
D a 



TD a ,M/3i 

i 

c-\ 



01 D a ,Mi3> 



Dg,M a - 



(2.47) 
(2.48) 
(2.49) 
(2.50) 
(2.51) 



Because the F-direction of Majorana particles is fixed to be equal to their A- 
direction, other vertices with Majorana particle having F-directions opposite 
to the A-direction do not appear in the G-lines. We also need the G-line rule 
for the fermion number violating chargino interactions with a conventional 
Dirac particle, whose lagrangian is given by 



Ccd = y c (C)T c , D V(D) + #(D)r AC * c (C), 



[2.52) 



where T D C = Iq^c ulo- The rules are 



C 
D 



D 

C 



—C 1 r c > j £ > , 



(2.53) 
(2.54) 



Table 2 and eq.(2.42)-(2.54) are what we need to compute the G-line ampli- 
tudes. 

As Denner et al. [16] have done, it is possible to formulate the fermion number 
clashing G-line amplitude without using the charge conjugation matrix C 
if one observes that the propagator with its F-direction opposite to the A- 
direction is given by 



-S T (-k) = -C^S^C, 



(2.55) 



k 



and that the right-hand side of (2.50) and (2.51) is expressed in terms of the 
charge conjugated vertex function 



(2.56) 
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as 



M/3 < • > D a r C M„ Da C, (2.57) 

!>,»— >•<— M a -C- 1 ^^, (2.58) 

In ref.[14] the equivalence of these two approaches has been demonstrated. 



3 Comments 

Several comments are in order. 

We have checked and debugged the possible coding errors in GRACE/SUSY by 
several methods. First, the interaction vertices coded in the system, which are 
based on the expression presented in [8], are compared with the result com- 
puted by the algebraic manipulation program REDUCE. As another check, we 
have used the gauge invariance of the processes which include photon, W ± , 
Z and gluon. In this check one has to make the widths of the heavy gauge 
bosons vanish. By changing the value of gauge parameters, we have confirmed 
at one phase space point that the amplitudes of all the possible 582,102 pro- 
cesses with up to 6 external legs remain unchanged within the error of 1 in 
10 15 or less in double-precision, and 1 in 10 30 or less in quadruple-precision 
calculation. Also further checks have been made by Boudjema and Belanger 
of LAPTH [17], who compared the GRACE/SUSY results of several SUSY pro- 
cesses with those of CompHEP [18], the automatic computation program for the 
high energy processes developed by a group of Moscow State University. 

Using GRACE/ SUSY we have computed a number of SUSY processes and some 
of the results have been published. As this system allows us to deal with any 
SUSY particle production process at various types of high energy colliders, 
there are analyses for e + e~ [5,19-22], ep [23-26], cy, 77 [27-30] and pp col- 
liders [31]. Note, moreover, that some works have been done with a version 
of GRACE/SUSY extended to the iZ-parity violating SUSY models [24-26,28- 
30,32]. 

GRACE/SUSY can treat as many final-state particles as possible, provided that 
one prepares by oneself the kinematics routine necessary for the phase space 
integration and the computer has an enough capacity to store the gigantic 
output files. The system is, however, equipped with its own kinematics library 
which has all the kinematics routines necessary for the analysis of up to 4-body 
final-state (see Appendix C). 

Users can modify the program for their own interest but at their own risk. In 
particular, if one would like to use the other set of input parameters for SUSY 
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lagrangian, for example, the trilinear coupling constants Af in place of Of, one 
should prepare a program in which the parameter values are converted from 
the user's set to that of GRACE/ SUSY given in Table 1, and call the program 
before the amplitude calculation starts. Note that in the system it is rrifAf 
but not Af itself that has the variable name (see Appendix B). Note also 
that in the sfermion sector, the unmixed physical state is obtained by setting 
either Of = or Of = it/ 2. In the former case, however, the identification 
corresponds to f\ = fi and f 2 = /r, while in the latter fi = f R and f 2 — ft,- 
Since < is assumed in the system, we have to set Of = (Of = n/2) 
when one considers the case in which there is no mixing in the sfermion sector 
and the left-handed sfermion is lighter (heavier) than the right-handed one. 

The program GRACE/ SUSY is designed strictly within the framework of the tree 
level. To be more realistic, however, one can arrange the value of the input 
parameters so that they are subject to the renormalization group equations 
(RGE). A program MUSE developed by Lafage [20,33] can compute the SUSY 
particle masses as solutions of the coupled RGE. One can combine such a 
program with GRACE/SUSY to invoke improved values for the input parameters. 
See also ref.[34]. 

It is also possible to extend GRACE/ SUSY so as to include new particles or new 
interactions such as the i?-parity violating interaction (see, for example, [26]). 

The algorithm of G-line works without any change for loop diagrams. Hence 
the extension of GRACE/SUSY for one-loop process is possible once one had 
appropriate loop evaluation routines. 

In its current version of GRACE/SUSY the quark-mixing matrix is taken to be 
the unit matrix. Neutrinos are massless and left-handed. One may introduce 
the quark mixing by modifying the model files in an appropriate way. 



4 How to use GRACE/SUSY 

In this section we show how to use the program GRACE/SUSY by adopting a 
SUSY process e + e~ — > iXiXi as an example. This process has 28 Feynman 
graphs at tree level, among which 4 graphs contain even number of fermion 
number violating vertices. 

4-1 Definition of physical process 

First the physical process to be computed must be defined, which is given 
in the file "in. pre". The content of the file in. pre for the process e + e~ — > 
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7Xi Xi 1S as follows; 

Model="mssm.mdl" ; 
Process ; 
ELWK=3; 
QCD=0; 

Init ial={electron , positron} ; 

Final ={photon, charginol , anti-charginol} ; 
Kinem="2302"; 
Pend; 



4.1.1 Model file 

In the first line the model file is defined. Two model files are prepared, one is 
mssm.mdl for the minimal SUSY standard model and the other is ms .mdl for 
the standard model. The latter will be useful when one wants to separate the 
SUSY contributions from those of non-SUSY in such processes where only the 
SM particles appear in the external legs. 



4-1-2 Order of perturbation 

The orders of perturbation are to be defined; ELWK for the electroweak inter- 
action and QCD for the QCD. 



4-1.3 Initial and final states 

In the lines for the initial and final states, the particle names of the process 
are to be given for each state according to the list given in Appendix A. One 
should be careful about the ordering of the particles (in particular, of massless 
gauge bosons) in each state. The program reads this input list and numbers 
all the particles from 1 to n for n external legs according to the order in which 
they appear in these two lines. These assigned numbers are in turn used to 
specify the integral variables in the phase space integration by the kinematic 
library. For example, for the 3-body final state including a single photon, the 
photon must appear first in the final state specification as our example shows. 
Otherwise the best choice of the integral variables will not be taken. See the 
next subsection and Appendix C for detail. In our example of e + e~ — > ^XiXi 
the input file in. pre defines the particle numbers 1~5 as e~ (pi), e + (p 2 ), 7 
(Ps), Xi (Pa) and Xi (Pb)- 
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4-1-4 Kinematics code number 

GRACE/ SUSY has a built-in kinematics library which contains a set of kinemat- 
ics routines for the phase space integration. In order to specify the kinematics 
routine to be used, one has to choose an appropriate code number, 2302 in 
our example. A list of code numbers is given in Appendix C. 



4-2 Graph generation 

The graph generation procedure starts by typing the command grc: 
grace /, grc 

All Feynman graphs for the process are generated and their information is 
saved in the file out . grf . 



4-3 Drawing graphs 

A Feynman graph drawer gracef ig is initiated by the command: 
grace / gracef ig 

The procedure gracef ig reads the graph information from the file out. grf 
and analyzes the structure of each graph, and then creates two windows. The 
output of our example process is shown in Fig.l. 

4-3.1 Drawing window of gracef ig 

On the drawing window the first m x n Feynman graphs are drawn, in the 
unitary gauge by default, where the number of rows m and the number of 
columns n are determined automatically by the size of screen and the total 
number of graphs. When the size of the window is changed, the number of 
columns n is conserved while the number of rows m is adjusted so as to fit the 
height of the window. 

4-3.2 Control window of gracef ig 

In the control window several buttons are set with which one can select through 
the pull-down menus several options which are listed in Table 3. 
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gracefig 




GRACEFIG : Cedetiby SJCaimtaec IHm&mttgsyc, KEK, JAPAN 



Co variant : 
Unitary 
Selected 

* mode : 
Draw Unitary 
1 -th out of 24 

* EPS -out : files 



Fig. 1. Feynman graphs for e + e — > 7Xi~Xi i n the MSSM. 

Particularly one can draw all the Feynman graphs in the covariant gauge by 
choosing an appropriate button on the menu. 



4-3.3 To draw selected graphs 

There are two modes, "Drawing Mode" and "Select Mode" . When gracefig is 
invoked, all Feynman graphs in the unitary gauge are displayed in the drawing 
mode. If one wants to see or print only some specific graphs, one can select 
them in the select mode by clicking their pictures in the drawing window. The 
background color of the selected graphs are darkened in this mode. After this 
selection only the selected graphs can be displayed in the drawing mode by 
choosing "Selected graphs" in the graph menu. 
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Table 3 



Buttons and sub-menus of gracef ig 



Button / Sub-menu 


Function 


Quit 


Exit from gracef ig. 




Next page 


Show next page. 




Previous page 


Show previous page. 




+5 page 


Show 5-th page ahead. 




-5 page 


Show 5-th page back. 




Scale up 


TV/Tl j 1 • C 11 mi 1 c 

Make the size of graph larger. Ihe number of j 
in a row: n—*n—l. 


graphs 


Scale down 


Make the size of graph smaller. The number of j 
in a row: n — > n + 1. 


graphs 


Option menu 


Pull down menu. 




/ Graph number 


On/Off display of graph number. 




/ Particle number 


On/Off display of particle number. 




/ Vertex number 


On/Off display of vertex number. 




/ Particle name 


On/Off display of particle name. 




/ Line number 


On/Off display of internal line number. 




Graph menu 


Pull down menu. 




/ Covariant Gauge 


Show graphs in covariant gauge. 




/ Unitary Gauge 


Show graphs in unitary gauge (Default). 




/ Selected Graphs 


Show only the selected graphs. 




Mode Selection 


Pull down menu. 




/ Drawing Mode 


Default mode to draw the graphs. 




/ Select Mode 


To select a graph, click on a graph you choose. 




EPS output 


Output the current page on a file grcf ig* . eps in eps 




format. 





4-3.4 EPS output and insertion of a picture in a text 

When "EPS output" button is clicked, the EPS file for the current page in the 
drawing window is created. The name of EPS file is "grcfignn.eps", where 
nn is the sequential number starting from 0, and one EPS file corresponds to 
one page of the graphs. 
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I BASES | | SPRING | 

mainbs mainsp 



bsinit 



userin 



BASES 



BSINFO ' 



BHPLOT ' 



BSWRIT ' 



setmdl 



setcpl 



sethptal 



gf init 



kinit 1 
i 

kfill ! 



f unc 



kinem < 



ampsum 



amptbl 



X 



amps 



X 



ampord 



Interface and 
CHANEL Library 



BSREAD 



SPRING 



' SPINFO 



' SHPLOT 



Fig. 2. Relation among the program components. 



4-4 Source code generation 



After the graph generation, by typing the command 
grace /, grcfort 

FORTRAN code and Makefile are generated, which are necessary to perform 
gauge independence check, the numerical integration and the event generation. 
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4-4-1 Generated source code 



Three kinds of program components are created by the system. They are used 
in the amplitude calculation, in the phase space integration by BASES and in 
the event generation by SPRING, respectively. 



The interrelation among the subprograms generated by GRACE/SUSY is de- 
picted in Fig.2, where those subprograms in the solid box are automatically 
generated by GRACE/ SUSY, while those in the dashed box are already contained 
in other program packages BASES/SPRING, kinematics library, interface pro- 
gram library to CHANEL and program package CHANEL. This figure omits a few 
minor components. The components vertically located in the middle of the fig- 
ure are commonly used in BASES and SPRING. The program specifications of 
BASES/SPRING, the interface and CHANEL are described in references [2-4,35]. 
A brief description of generated subprograms is given below: 

(1) Main programs. 



mainbs 


( main ) 


is the main program for the integration. 


mainsp 


( main ) 


is the main program for the event generation. 


gauge 


( main ) 


is the main program for the check of gauge 






independence. 



(2) A set of program components for the integration by BASES and for the event 
generation by SPRING. 



userin 


( subroutine ) 


controls initialization. 


setmdl 


( subroutine ) 


defines MSSM parameters, masses and decay 
widths of particles. 


sethptal 


( subroutine ) 


defines the parameters of the extended Higgs 
potential. 


setcpl 


( subroutine ) 


defines coupling constants. 


gf init 


( subroutine ) 


initializes particle polarization and selects Feyn- 
man graphs (amplitudes). 


kinit 


( subroutine ) 


initializes integration parameters, histograms, 
kinematics and user's parameters. 


f unc 


( function ) 


calculates the numerical values of differential 
cross section. 


kinem 


( subroutine ) 


derives particle four momenta from the integra- 
tion variables. 


kf ill 


( subroutine ) 


fills values in histograms and scatter-plots. 
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amptbl 


( subroutine ) 


calls amnnnn to calculate amplitudes. 


amp sum 


( subroutine ) 


sums the square of the matrix elements 






cwpr trip nplipitv ^t^itp^ 


amps 


( subroutine ) 


calculates amplitudes. 


ampord 


( subroutine ) 


arranges amplitudes. 


incll .h 


( include file ) 


defines the common variables for masses, am- 






"nlitnrlp farilpG c>~t~f 
piltULIC LdUlCo OIL. 


incl2.h 


( include file ) 


defines the work space for amptbl. 


inclk.h 


( include file ) 


transfers values of masses and a few constants. 


userout 


( subroutine ) 


prints the amplitude summary table. Used only 






by gauge. 



The subprograms userin and f unc are used both in the numerical integration 
by BASES and in the event generation by SPRING. In userin the subroutine 
kinit is called to initialize the kinematics routine and other parameters. The 
subroutines setmdl, setcpl and gf init are also called for the initialization 
of physical parameters and amplitude calculation. 



The function subprogram func is used for calculating the numerical value 
of differential cross section, where the subroutine kinem is called for calcu- 
lating four-momentum vectors of external particles, and amptbl and ampsum 
are called for the amplitude calculation. The detailed specifications of these 
program components are described in refs. [2,35]. 

4-4-2 Setting physical parameters for numerical calculation 

Prior to running the program, one has to set the values of the input physical 
parameters by editing appropriate FORTRAN source codes. The parameters and 
corresponding subroutines to be edited are categorized as follows. 

kinematical parameters kinit, kf ill 
MSSM parameters setmdl, sethptal 
polarization and graph selection gf init 

(1) Initialization routine kinit 

The values of the following parameters can be changed by editing kinit. 
(a) Modification of physical parameters 

Center of mass energy in GeV. ^/s = lOOOGeV is a default setting, 
which corresponds to a description w = 1000. dO in kinit. 
Various kinematical cuts etc. 
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Option flags in built-in kinematics routines ([36]). 

(b) Integration parameters which control the action of BASES can be altered in 
kinit. The following parameters can be changed from their default values 
according to the required integration accuracy: 

ncall The number of sampling points per iteration 

accl The expected accuracy for the grid optimization step 

itmxl The maximum iteration number for the grid 
optimization 

acc2 The expected accuracy for the integration step 

itmx2 The maximum iteration number for the integration step 
However, the following parameters are usually not to be changed: 
ndim The dimensions of integral 

nwild The number of wild variables 

xl(i) ,xu(z) The lower and upper bounds of integration variable x(i) 

ig(0 The grid optimization flag for i-th variable 

(c) Initialization of histograms and scatter-plots 

If the user wants histograms for some physical variables, the initializa- 
tion routines xhinit and dhinit should be called for each histogram and 
scatter-plot, respectively. 

The system automatically generates source code to call these routines 
so as to construct histograms with respect to each integration variable, 
x (i) , i=l ~ ndim and the energies and angles of final particles. If the user 
does not need these default histograms, one may delete the corresponding 
FORTRAN statements in kinit. When the user needs more histograms for 
some distributions, one should add lines to call xhinit or dhinit. These 
routines are to be written in the following format 
call xhinit ( id#, 



where 1 < id# < 99. 
(2) Filling histograms and scatter-plots in kf ill 

In the generated kfill by GRACE/SUSY, histograms for all integration 
variables and scatter-plots for all combinations of them are filled by calling 



lower_limit , upper_limit, # of bins, 
'Title of this histogram' ) 



and 

call dhinit ( id#, 



x-lower_limit , 
y-lower_limit , 
'Title of this 



x-upper_limit , # of x bins, 
y-upper_limit , # of y bins, 
scatter plot' ) 
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xhf ill and dhf ill, respectively. If one changes the initialization of his- 
tograms and scatter-plots in kinit, their filling parts should be changed 
accordingly. These routines are to be called as 

call xhf ill ( id#, v, func ) 

call dhf ill ( id#, vx, vy, func ) 

(3) Subprogram setmdl 

The MSSM parameters, masses and widths of particles are defined in 
setmdl subroutine. The variable names of all the parameters are listed in 
Appendix B. As an example we show how to set the MSSM parameters in 
the case (//, M 2 , tan/5, m u ~ e ) = (-300GeV, 200GeV, 10, 150GeV). 

tanbe = 10. dO 

xmu = -300 . dO 

xm2 = 200 . dO 

amsn(l) = 150. dO 

(4) Subprogram sethptal 

In this subprogram one has to specify the Higgs potential to be used 
in the calculation by assigning a value 1 or 2 to the variable ihiggs. See 
section 2.4 for detail. 

■ ihiggs = 1: mode 1; MSSM Higgs potential 

• ihiggs = 2: mode 2; general extended Higgs potential (default) 

The default value is ihiggs=2, in which case one has to give values of the 

additional seven input parameters, for example, 

amhl = 115. dO 

amh2 = 300. dO 

amhc = 350. dO 

al = -0.37d0 

c4 = O.dO 

c9 = O.dO 

clO = O.dO 

(5) Subprogram gf init 

(a) Selection of graphs. 

The contribution from a graph is calculated when its flag is set to be 
jselg = 1 in gfinit. For example, if one wants to know the contribu- 
tions only from the sneutrino exchange (graphs 19, 20, 21 and 28 in Fig.l), 
one should write 

do 10 nl = 1, ngraph 
jselg(nl) = 
10 continue 

jselg(19) = 1 
jselg(20) = 1 
jselg(21) = 1 
jselg(28) = 1 
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The selection of graphs should respect gauge invariance, otherwise the 
cross section is physically not sensible. It should be noted that this speci- 
fication of graphs is independent from the graph selection in the Feynman 
graph drawer gracef ig. 
(b) Set of polarization. 

As the helicity amplitude technique is used in GRACE/ SUSY, the beam 
polarization can be specified. The polarization sum is performed by set- 
ting the flags jhs(i) and jhe(i) for particle i having spin Sj, where 
< jhs(i) < jhe(i) < 2sj. In default the unpolarized cross section is 
assumed, which corresponds to setting the flags for the initial e~ (particle 
1) and e + (particle 2) as 

* 1: initial electron mass=amlp(l) 



jhs(l) = 
jhe(l) = 1 



* 2: initial positron mass=amlp(l) 



jhs(2) = 
jhe(2) = 1 



When we consider the initial e beam polarization, for example, they 
should be changed to 

jhs(l) = 

jhe(l) = 

for helicity=-l and 
jhs(l) = 1 
jhe(l) = 1 

for helicity=+l. 

Similarly, one can choose polarization vectors of external gauge bosons by 
putting jhs (i) , jhe(i) = 0, 1 (transverse), 2(longitudinal: only for mas- 
sive gauge bosons). In the default the linearly polarized states are realized 
as combinations of the transverse polarizations. When one wants to take 
the circular polarization for a photon, an option switch jcpol should be 
changed from to 1. 
(c) Gauge parameters. 

Without a specification the unitary gauge is used in the integration over 
the phase space. However, one can select covariant gauge (i?g-gauge) with 
the gauge parameters having arbitrary values. 

The selection of gauge is controlled by integer variables igauwb, igauzb, 
igauab and igaugl which correspond to gauge bosons W ± , Z, 7 and gluon, 
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respectively. When these parameters are set to be 0, unitary gauge is se- 
lected. Otherwise they should be igauwb = 1, igauzb = 2, igauab = 3 
and igaugl = 4. The value of each gauge parameter should be given by 
agauge (igauxx) , where xx is one of wb, zb, ab or gl. 

In the file gf init . f one will find the following lines which select the 
unitary gauge. 

* Gauge parametes (default is unitary gauge) 

igauab = 
igauwb = 
igauzb = 
igaugl = 

agauge (0) = 1.0d20 
In a general covariant gauge these lines should be changed as follows; 

* Gauge parametes (covariant gauge) 

igauwb = 1 
igauzb = 2 
igauab = 3 
igaugl = 4 

agauge (igauwb) = ... 

agauge (igauzb) = ... 

agauge (igauab) = ... 

agauge (igaugl) = . . . 

4-5 Compilation and gauge independence check 

As described in a previous subsection, Makefile is also generated by the 
command grcf ort. The command make creates the executable files gauge, 
integ and spring. 

Before calculating the cross section a preliminary test of the program source 
should be done. In GRACE/SUSY system gauge independence test is prepared for 
this purpose. The numerical values of the differential cross section at one sam- 
ple point in the phase space are calculated both in the covariant and unitary 
gauge, and the consistency between these two values is examined. By exe- 
cuting gauge, the contributions to the cross section at one phase space point 
in unitary gauge (ansl) and in covariant gauge (ans2), and their difference 
(ansl/ans2-l) are shown as the following example demonstrates. 



ansl = 0. 1255918104665936E+01 
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ans2 = 0.1255918104665936E+01 



ansl/ans2 - 1 = 2 . 22044605E-16 



The relative difference should be around 10 16 in the double-precision com- 
putation and around 1CT 30 in the quadruple-precision. 

It should be noted that even though this gauge independence check shows con- 
sistency between the two gauges, it does not necessarily guarantee a complete 
gauge independence, since it tests only at a specific point in the phase space. 
It is recommended to repeat a similar check at several points in phase space 
by changing the variable x(i) in gauge. f from their default value of 0.45 to 
other value lying in the range < x(i) < 1. 

4-6 Integration 

After a successful test of gauge independence, one can proceed to the numer- 
ical integration by BASES, whose details can be found in ref.[2]. By the make 
command one has already created the executable files integ for the integration 
and spring for event generation. 

For integration the command integ is used. 

grace /, integ 

The integration parameters and the result of the integration at each iteration 
are shown on the screen, and at the same time they are automatically stored 
in the file bases. result together with the histogram outputs. If one uses a 
built-in kinematics routine, the cross section is given in unit of pico-barn (pb). 

The numerical integration is performed by two successive steps, the first is 
the grid optimization step and the second the integration step. In the former 
the gridding of the phase space in hypercubes is optimized for the integrand, 
while in the latter an accurate estimate of the integral is evaluated with frozen 
grids determined in the first step. 

Before termination of the integration procedure, BASES creates a binary file 
bases. data and writes the probability information on it, which is used for 
the event generation. 

The system will issue WARNING messages when the convergence is not well es- 
tablished. It is always recommended to look at the integration result carefully, 
particularly over the convergence behavior both for the grid optimization and 
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integration step. When the convergence accuracy fluctuates largely from one 
iteration to another, and/or, when in some cases it suddenly jumps up to a 
value far from the preceding iterations, the resultant estimate of integral is 
not reliable. There are two possible origins of such a behavior; one is due to 
too small number of sampling points and the other, more serious, is due to an 
unsuitable choice of the integration variables for the integrand. In the latter 
case one should use a better set of the integral variables. 

4-7 Event generation 

The event generation starts with the command 
grace"/, spring 

Then SPRING reads the probability information from the binary file bases . data 
and asks the number of events with the following prompt: 

Number of events ? 

The user must type the number of events to be generated. 

If the sample point does not meet the generation criteria, such a point is 
discarded as a "failed generation". The event generation will run until a given 
number of events are generated or the number of failed generations exceeds 
its given maximum mxtry, whose default value is 50. 

In order to estimate the computing time for the event generation, one can refer 
to the expected generation time given in the information of BASES output, 
where the generation efficiency of 70% is assumed. 

4-7.1 Analysis of events 

In the event generation four-momentum vectors of the final particles are au- 
tomatically stored in vec(i,n) in mainsp, where n stands for the particle 
number and i = 1 ~ 4 correspond to the three- momentum (1, 2, 3) and the 
energy (4) of the particle. 

To create histograms of some variables of the generated events, one has to ini- 
tialize them before the generation loop and fill them after each call of spring. 
This kind of histograms is called as an " additional histogram" in contrast to 
the " original histogram" which is filled both in the integration and event gen- 
eration steps. In the original histograms one can see how the generated events 
reproduce those distributions produced by the integration, while in the addi- 
tional histogram only the frequency distribution can be seen as in the usual 
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histogram. An example for saving the four-momentum vectors and creating 
additional histograms is given below. 

* Open the output file for the events 



* Initialize the additional histograms by 

call xhinit( id, ... ). 

do 100 nevnt = 1, mxevnt 

call spring ( func, mxtry ) 

* 

* analyze the event and 

* fill the additional histograms by 

call xhfilK id, ) . 

* 

do 90 k = 1 , nextrn 

write(6,*) (vec(j ,k) , j=l ,4) 
90 continue 
100 continue 

4.7.2 Output from SPRING 

The output from SPRING is written on the file named spring. result, which 
consists of the general information, original and additional histograms, scatter- 
plots, and the number-of-trials distribution. In the last distribution one can 
see the generation efficiency. It is recommended to check if this distribution 
has a sharp peak in the first bin. If this is not the case, the integration should 
be carried out again with more sampling points or with more iterations. 
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A List of particles in MSSM 



The particles appearing in MSSM and their variable names are listed in two 
tables below. These variable names are to be used in specifying the process in 
in. pre. 



particle 


variable name 


photon 


photon 


w + {w~) 


W-plus (W-minus) 


z 


Z 


gluon 


gluon 




nu-e (nu-e-bar) 


e-(e+) 


electron (positron) 




nu-mu (nu-mu-bar) 




muon (anti-muon) 




nu-tau (nu-tau-bar) 


r (r+) 


tau (anti-tau) 


u(u) 


u (u-bar) 


d{d) 


d (d-bar) 


c(c) 


c (c-bar) 


8(8) 


s (s-bar) 


W 


t (t-bar) 


6(5) 


b (b-bar) 


h° 


Higgsl 


H° 


Higgs2 


A 


Higgs3 


H+{H-) 


Higgs-plus (Higgs-minus) 
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particle 


variable name 


XiiXi) 


charginol (anti-charginol) 


xi(x2 ) 


chargino2 (anti-chargino2) 


Xi 


neutralinol 


x\ 


neutralino2 


x% 


neutralino3 


x\ 


neutral ino4 




snu-e (anti-snu-e) 




snu-mu (anti-snu-mu) 


z7 r (z7 T ) 


snu-tau (anti-snu-tau) 


e~r(el) 


selectronl (anti-selectronl) 


e 2~( e 2 ) 


selectron2 (anti-selectron2) 




smuonl (anti-smuonl) 


^2(^2) 


smuon2 (anti-smuon2) 


t~i(t?) 


staul (anti-staul) 


r~2~{r~t) 


stau2 (anti-stau2) 


U\{ui) 


sul (anti-sul) 


4(4) 


su2 (anti-su2) 




sdl (anti-sdl) 


4(4) 


sd2 (anti-sd2) 


c~i(ci) 


scl (anti-scl) 




sc2 (anti-sc2) 


J- V 1 / 


ssl (anti-ssl) 


59(52) 


ss2 (anti-ss2) 




stl (anti-stl) 




st2 (anti-st2) 


&i(&i) 


sbl (anti-sbl) 


& 2 (5) 


sb2 (anti-sb2) 
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B Name List of GRACE/SUSY Parameters 



The variable names used in GRACE/SUSY are summarized. All mass variables, 
widths and the dimensional coupling constant \i must be given in unit of GeV. 
Those variables whose magnitude can be specified by the user in setmdl.f 
and sethptal.f are indicated by an asterisk "*" at the end of the row. The 
value of other variables should not be changed by the user. Such a 
modification may cause internal inconsistencies. 



gauge bosons 



category 


notation 


variable name in GRACE/ SUSY 




gauge boson masses 


M z , M w 


amz, amw 


* 


gauge boson widths 




agz, agw 


* 


coupling 




alpha, alphas 


* 




e, g s 


ge, gs 






g, g' 


gg, g z 




mixing angle 


cos#w, sin6*vK, tan#vK 


cw, sw, tw 





Ordinary fermions 



fermion mass 


m u , m c , m t 


amuq(l), amuq(2), amuq(3) 


* 




m d , m s , m b 


amdq(l), amdq(2), amdq(3) 


* 




m e , m M , m T 


amlp(l), amlp(2), amlp(3) 


* 


fermion widths 




aguq(l), aguq(2), aguq(3) 


* 






agdq(l), agdq(2), agdq(3) 


* 




r r r 

1 ei ill t 


aglp(l) , aglp(2), aglp(3) 


* 



Of course one should set aguq(l) = agdq(l) = aglp(l) = 0. 



Higgs bosons 
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SUSY parameters 


(j, 


xmu 


* 


Higgs masses 


M A o 


amh3 


* 




M h ,M H , M H ± 


amhl, amh2, amhc 


* 


Higgs widths 


r\, t h , r A , r H ± 


aghl, agh2, agh3, aghc 


* 


mixing angle 


a, cos a, sin a 


al, ca, sa 


* 




cos f3, sin /3,tan (3, cot f3 


cb, sb, tanbe, cotbe 


* 



Charging 



SUSY parameters 


A*, M 2 


xmu, xm2 


* 


chargino masses 


m xt> m xf> 


amsw(l), amsw(2) 




chargino widths 


T xt> T xt> 


agsw(l), agsw(2) 


* 


chargino mixing angles 


sin0 L , cos0 L , e L 
sm<f) R , cos(p R 


sphl, cphl, epsl 
sphr, cphr 





Neutralino 



SUSY parameters 


fi, Mi, M 2 


xmu, xml, xm2 


* 


neutralino masses 




amnl(i) 




neutralino widths 




agnl(i) 


* 


neutralino mixing matrix 


(0 N )ij, rji 


ogmn(i,j), etan(i) (complex) 





Gluino 



gluino mass 


m~ 9 = M 3 


amgl = xm3 


* 


gluino width 




agg 1 


* 



Sfermion 
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sfermion masses 


m,-;., Trip. , rrif 


amsu(l,i), amsu(2,i), amsu(3,i) 


* 




my., m Si , mi. 


amsd(l,i), amsd(2,i), amsd(3,i) 


* 




m De , m^, mo T 


amsn(l), amsn(2), amsn(3) 


* 




m e -,, rrifc, m fi 


amsl(l,i), amsl(2,i), amsl(3,i) 


* 


sfermion widths 


V- r~ r~ 

1-11.. ±. ji . . J.J.. 


agsu(l,i), agsu(2,i), agsu(3,i) 


* 




r~ r~ r~ 


agsd(l,i), agsd(2,i), agsd(3,i) 


* 




r~ r~ r~ 


agsn(l), agsn(2), agsn(3) 


* 




V- V- V- 


agsl(l,i), agsl(2,i), agsl(3,i) 


* 


SUSY-breaking 


m u A u , m c A c , m t A t 


xauq(l), xauq(2), xauq(3) 




parameters 


m d A d , m s A s , m h A b 


xadq(l), xadq(2), xadq(3) 






m e A e , m^A^, m T A T 


xalp(l), xalp(2), xalp(3) 




sfermion mixing 


sin^u, sin^c, angles sin^ 


shuq(l), shuq(2), shuq(3) 


* 




COS#«, COS 9c, COS#£ 


chuq(l), chuq(2), chuq(3) 


* 




sin #7, sin^s, shift- 


shdq(l), shdq(2), shdq(3) 


* 




COS 9 J, cos cos 6^ 


chdq(l), chdq(2), chdq(3) 


* 




sin 9g, sin Op,, shift- 


shlp(l), ship (2), ship (3) 


* 




COsftejCOsft^, COS#f 


chlp(l), chip (2), chip (3) 


* 



C Kinematics Library 

C.l Cross section formula 

The volume element for each of final particles is given by 

d 3 p 



2£(2tt) 3 ' 

where p is the three momentum of the particle and E is its energy. 
In a collision of two particles labeled by 1 and 2, the flux is given by 

Pi P2 

flux = v rel 2E 1 2E 2 , v rel = — - — 
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(C.l) 



(C.2) 



with v re i being the relative velocity of particles 1 and 2. The total cross section 
is then written as 

* = ik n / +P2-Eft)EE Pi,i', (c.3) 

J 7 3 h f hi 

where Tjj is the helicity amplitude for the process i — > /. Initial helicities are 
averaged and the final ones are summed. This formula is used for the cross 
section in the GRACE/SUSY system. The cross section is given in the unit of 
pico-barn(pb) in the kinematics library. 

For the decay of a particle with momentum p, the total decay rate is given by 

The independent integration variables expressed in terms of momenta and 
angles must be mapped onto those of the integration package BASES which are 
normalized as < x(i) < 1, where i runs from 1 to ndim, the dimension of 
the integration, 

One should remember that the statistical factor for the identical 
particles is not generated automatically and should be supplied by 
hand. This is because this factor depends on the observables one is interested 
in. For the total cross section one should divide the result obtained by BASES 
by the factorial of the number of identical particles. 



C.2 Kinematics library 

GRACE/SUSY assumes the following two types of processes: 

Pi + P2 — > + Pa H (2-body scattering), 

Pi —> P2 + P3 + Pa + • • ' (decay of a particle). 

It should be noted that the particle assignment follows the order in which the 
particles are originally defined in in. pre. 

The Lorentz invariant phase space element for a final n-body (A — > 1 + 2 + 
3 + \-n) is defined by 
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df n = (2tt) 4 <5 4 (V> - T,p\ TT 7 ^ = 7 I 7 dF n , (C.5) 

^„=5 4 (Ep-Ep)f[^- (C6) 

\ in out / out j 



The following chain rule is useful(0 < k < n — 1): 



rfr„(A ->i + 2 + 3 + ---) 

dO 2 

= dT k+1 (A^l + --- + k + q)-^- 

In 

xdF n _ k (q^ (fc + l) + (fc + 2) + •••), (C7) 



where g 2 = Q 2 



When p a is in the center-of-mass system, i.e., p a = 0, and if p a = Pb + p c we 
introduce the 2-body phase space 



dT CM (a; be) = dY 2 = 8\p a - p b - Pc)^^, (C.8) 



and write it by angular variables: 



dT C M{a; be) = ^-^ — g^cm(«; 6c) = ^ Q ' &C dcos9 b ^ bc )d(f) b ^ bc ). (C.9) 



Here /5 is given by 



2P 

/3(a;6c) = — (CIO) 



v (E'a + m fe + m c )(E a - m b - m c )(E a + m b - m c )(E a -m b + m c ). 



The subscript (6c) of the angles indicates that 9 b , <f> b are defined in the center- 
of-mass frame. Angles in the laboratory frame have no subscript. 
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Kinematics routines prepared in the library are summarized in Table C.2. 



code number 


contents 


1201 


1-body — > 2 body decay 


1301 


1-body — > 3 body decay 

Sequential decay 1^2 + (3 + 4) ^2 + 3 + 4 can be treated. 


2201 


2-body — > 2 body in CM frame 

t- and -u-channel singularities can be treated. 


2301 


2-body — > 3 body in CM frame , 

Sequential decay type 1 + 2-^3 + (4 + 5) ^3 + 4 + 5. 
Resonance on particles 4 and 5 can be treated. 


2302 


2-body — > 3 body in CM frame , 
Radiative processes 1 + 2 — > 3(7) + 4 + 5, 
both initial and final radiation can be treated. 


2303 


2-body — > 3 body in CM frame , 
Double-radiative processes 1 + 2 — > 3(7) + 4(7) + 5 


2304 


2-body — > 3 body in CM frame , 

Three photon processes 1 + 2 — * 3(7) + 4(7) + 5(7) 


2401 


2-body — > 4 body in CM frame, a pair of sequential 
decay type 1 + 2^ (3 + 4) + (5 + 6) ^3 + 4 + 5 + 6 
t-channel singularity can be treated. 


2402 


2-body — > 4 body in CM frame, 

'fusion' type 1 + 2 -> (3 + A) + (4 + B); A + B -> 5 + 6 



Table C.2 Kinematics library prepared in GRACE/SUSY system. 
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